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ABSTRACT 

Steady accretion of a radiating gas onto a central mass point is described and compared to classic 
Bondi accretion. Radiation losses are essential for accretion flows to be observed. Unlike Bondi 
flows, radiating Bondi flows pass through a sonic point at a finite radius and become supersonic near 
the center. The morphology of all radiating flows is described by a single dimensionless parameter 
proportional to M/MT S where T s is the gas temperature at the sonic point. In radiating Bondi flows 
the relationship between the mass accretion rate and central mass, M oc M p with p ~ 1, differs 
significantly from the quadratic dependence in classical Bondi flows, M oc M 2 . Mass accretion rates 
onto galaxy or cluster-centered black holes estimated from traditional and radiating Bondi flows are 
significantly different. In radiating Bondi flows the gas temperature increases at large radii, as in the 
cores of many galaxy groups and clusters, allowing radiating Bondi flows to merge naturally with gas 
arriving from their cluster environments. Some radiating flows cool completely before reaching the 
center of the flow, and this also occurs in cooling site flows, in which there is no central gravitating 
mass. 

Subject headings: hydrodynamics, galaxies: cooling flows, galaxies: clusters 



1. INTRODUCTION 

Bondi flow describes the steady spherical inflow of an 
adiabatic gas toward an accreting mass point. As gas 
flows radially toward the gravitating mass, its density in- 
creases due to tidal compression by the gravity field and 
consequently its temperature also increases. Gas at rest 
at infinity with uniform temperature and density acceler- 
ates steadily as it flows inward toward the central mass. 
In Bondi flow the mass of flowing gas is assumed to be 
small so its self-gravity can be neglected and the central 
mass can be regarded as constant during relevant flow 
times (Bondi 1952). Conventional non-radiating Bondi 
flows are often used in computational studies of massive 
black hole growth during galaxy formation or to detect 
errors in hydrocodes. Bondi flows are also often used 
to estimate black hole accretion rates from X-ray obser- 
vations of gas near central black holes in galaxies and 
galaxy clusters. However, Bondi flows do not radiate. 

For more realistic astrophysical applications, such as 
the flow of hot intracluster gas toward a cluster-centered 
black hole, the gas must radiate as it flows inward. Our 
purpose here is to illustrate that steady spherical accre- 
tion flow solutions with radiative losses differ fundamen- 
tally from the traditional Bondi solution. A large family 
of morphologically varied radiating solutions are possi- 
ble, all characterized by a single dimensionless parame- 
ter. But the most important difference is that the often- 
assumed quadratic Bondi relation between the mass ac- 
cretion rate and the central mass, M oc M 2 , no Ion ger 
holds. Instead, radiating flows require M oc M p with 
p ~ 1. Unlike classical Bondi flow, the gas temperature 
in radiating solutions always increases at large radii, al- 
lowing more natural fits to the hot gas cores in galaxy 
groups and clusters that also have positive temperature 
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gradients. Radiating Bondi flows can also accommodate 
larger accretion rates M due to gas flowing in from an 
extended cluster environment. 

To be useful in cosmological evolutionary computa- 
tions and in interpreting accretion rates from X-ray ob- 
servations, it is useful to keep our solutions as simple 
and general as possible. Consequently, our radiating 
steady solutions do not include spatially distributed mass 
of stars or dark matter (e.g. Yahil & Ostriker 1973; 
Cowie & Binney 1977; Fabian & Nulsen 1977; Mathews 
& Bregman 1978; Quataert & Narayan 2000; Guo, Oh, 
& Ruszkowski 2008), feedback energy, or rotation (Park 
2009; Brighenti et al. 2009; Narayan & Fabian 2011). 

When radiative losses are included, an initially hot gas 
can sometimes cool to zero temperature before it reaches 
the concentrated mass at the center, particularly if the 
central mass is sufficiently small or if there is no central 
gravitating mass at all. This latter type of flow toward 
a cooling site in the absence of gravity experiences no 
tidal compression so the gas temperature decreases con- 
tinuously toward the origin, becoming supersonic before 
cooling to zero temperature. Cooling site flows may de- 
scribe the appearance and evolution of cool, non-central 
low-entropy gas often observed in galaxy clusters far from 
the gravitating cluster core. 

2. FLOW EQUATIONS 

Consider the equations for steady spherical Bondi ac- 
cretion toward a point mass M but also include optically 
thin radiative losses: 

M = puiirr 2 (1) 
du dP GM 
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The mass flow M is negative for inflowing gas, u < 0. 
The specific thermal energy e = (P/p)/(7 — 1) = 3P/2p 
with 7 = 5/3. For simplicity we assume that the opti- 
cally thin cooling coefficient A (erg cm 3 s _1 ) is constant 
throughout the flow. While the temperature dependence 
of A is well-known (e.g. Sutherland & Dopita 1993), 
the quadratic density coefficient in the cooling term con- 
tributes much more to the non-linear cooling process 
than the comparatively slower variation of A with tem- 
perature. Moreover, a constant A is easy to implement 
in checking gasdynamical hydrocodes. 

The three flow equations can be reduced to two by 
eliminating the density p using equation (1): 
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These equations become singular at the sonic radius 
r = r s where the flow speed equals the local adiabatic 
sound speed, 

u s = c s = {59J3) 1 ' 2 . (8) 

The sonic radius is found by insisting that both square 
brackets in the equations above vanish when u = 
{56/3) 1 ' 2 
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Use of L 'Hospital's rule on the right side of either equa- 
tion (4) or (5) at the sonic point singularity gives the fol- 
lowing logarithmic slopes at the sonic radius (subscript 
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Beginning at the sonic point, equations (4) and (5) can 
be integrated inward and outward using the sonic slopes 
in equations (10) and (11). The sign option in equation 
(10) can be chosen for inflow or outflow. We restrict our 
discussion here to inflow solutions (minus sign). 

The structure of solutions of equations (4) and (5) is 
more apparent when transformed to dimensionlcss vari- 
ables 
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in which the differential equations become 
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It is remarkable that the morphology of all flow solu- 
tions ?7(£) and r(£) depends on only one dimensionless 
parameter 
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where p. — 0.61. Since M < 0, a must also be negative 
or zero and we consider values in the range — oo < a < 0. 

When radiative cooling is small, A — > 0, both a and 
A also become small and equations (4) and (5) describe 
standard adiabatic Bondi flow. When 7 = 5/3, which we 
consider here, the traditional non-radiating Bondi accre- 
tion solution is unique, entirely subsonic and singular 
since the sonic point lies exactly at the origin, r s = 0. 
Since our dimensionless variable £ = r/r s is inappropri- 
ate when r s = 0, we do not discuss the transition to 
standard adiabatic flow in detail here. When the central 
mass is small, a — > — 00, solutions approach mass- free 
flow toward a cooling site with M = 0. For sufficiently 
large |a| inflowing gas can cool completely (9 — > 0) before 
reaching the origin. 

3. RADIATING BONDI FLOWS 

The top and central panels of Figure 1 show eight su- 
perposed solutions of equations (14) and (15), ry(£) and 
t(£), each corresponding to a different value of the pa- 
rameter |a| = 10~ 6 , 10~ 2 , 0.1, 1, 1.5, 3, 10 and 100. The 
dotted extensions toward the left for each of these so- 
lutions describe supersonic flow within the sonic radius 
£ < 1. The dimensionless temperature r = 6/6 s , shown 
in the central panel of Figure 1, has the largest variation 
with \a\. For \a\ > 3 all plotted flows cool monotoni- 
cally to zero temperature at a finite radius ^ coo i not far 
within the sonic radius. At zero temperature the cooled 
gas would go into free fall toward the central mass, but 
this is not explicitly calculated. 

The parameter a is related to ratios between the flow, 
free fall and cooling times 

tff = Z[ a i/2 t Siow = f H iff = 7T \a\ 
4 t coo i 3 a tfcooi 3 a 1 / 2 



where a = 1 + (2/5) |a| and the times 

r s it r s °> £ s m p 

t flow = — tff = - (2GM)1/2 t coo l = -^- 

are evaluated at the sonic radius r s . 

The dimensionless velocity variable i] (£) = u 2 /u 2 de- 
creases monotonically with radius for all radiatively cool- 
ing solutions. For \a\ < 10~ 6 the top panel in Figure 1 
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shows that the velocity variable r](£) is nearly a single 
power law but becomes a double power law for \a\ > 0.1, 
broken at the sonic radius £ = 1. 

Of particular interest are radiating flows with \a\ > 3 in 
which the gas temperature drops to zero before it reaches 
the flow center. These solutions are expected when the 
central mass M (and/or the sonic point temperature T s ) 
becomes sufficiently small or when M is unusually large. 
When \a\ > 1, the overlapping dotted lines for 7j(£) in 
Figure 1 obscure the different leftward terminations of 
r/(£) which do not extend within the radius £ coo z where 
r — > in the bottom panel, with £ coo ; increasing with 
\a\. 

A temperature minimum at £ = £ m i„ characteristic of 
radiating Bondi flows having smaller \a\ appears in Fig- 
ure 1 for the first time at log£ m i„ w — 1 in the supersonic 
flow for \a\ = 1.5. As \a\ decreases further, the location 
£,min of the minimum in r(£) moves progressively toward 
larger radii: log£ m i„ « —0.1, 0.6, and 1.4 for \a\ = 1, 0.1 
and 0.01 respectively. For \a\ < 10~ 6 the temperature 
minimum has moved to log £ > 2 and no longer appears 
in Figure 1. For these very small values of \a\ the profiles 
for both r}(£) and r(£) converge to power laws in £ < 100 
that appear as single straight lines in the region plotted 
in Figure 1. 

Using equation (1), the corresponding dimensionless 
gas density variable 
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is normalized with the gas density at the sonic radius p s . 
The bottom panel in Figure 1 shows that decreases 
monotonically with radius for all solutions. 

The flow profiles in Figure 1 for velocity, temperature 
and density in dimensionless variables have morphologies 
that depend only on the single dimensionless parameter 
\a\. However, when these results are converted into phys- 
ical units, (two of) the fundamental dimensional flow pa- 
rameters in equation (16) - M, M, T s and a - become 
individually important, particularly for the normalizing 
factors r s and p s . The sonic radius 
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depends on GM/9 S as well as a and the density at the 
sonic radius 



M 1 / 3 



1/2 



(20) 



has its own unique dependence on flow parameters. 

All flows with \a\ ^ 3 have a central thermal peak. It 
is of interest to consider the total energy emitted in the 
central thermal peak £ < ^ m in 
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The asymptotic behavior of flow variables in the ther- 
mal peak as £ — > is: r\ cx £ _1 , 5 oc £~ 3 / 2 and 
r oc Consequently, the integral for the total peak 

luminosity, oc Jd£/£, diverges logarithmically at the ori- 
gin. (If A had a bremsstrahlung temperature dependence 
A oc T 1 / 2 , the divergence would be even stronger.) Evi- 
dently this divergence is related to the concept of a point 
mass which can increase the internal energy and density 
of radially inflowing gas without limit. Obviously rota- 
tion, magnetic or relativistic effects must eventually alter 
the nature of idealized radiating Bondi flows within some 
radius, removing the apparent divergence of the luminos- 
ity L from the thermal peak. 

In radiating Bondi flows having central thermal peaks 
the temperature in the peak increases as r oc but 
the dimensionless entropy factor, a = t/S 2 / 3 oc £°, is 
constant and the radiative cooling time is longer than 
the flow time. Gas flows supersonically to the center 
along an adiabat. But this adiabat differs from that of 
classical adiabatic Bondi accretion for 7 = 5/3 in which 
gas flows subsonically toward the central point mass, be- 
coming sonic as r — > r s = 0. In our radiating flows 
the rapidly converging and accelerating central flow be- 
comes effectively adiabatic for all values of a because 
of the short flow time, not because it stops radiating. 
Moreover, this central pseudo-adiabatic flow in r < r s 
is entirely supersonic, unlike the subsonic central flow in 
classical Bondi flow. 

The adiabatic character of asymptotic flow in central 
thermal peaks can be understood by comparing the two 
terms on the right in equation (3) that describe the im- 
portance of compression heating and radiation losses re- 
spectively. The ratio of thermal energy change due to 
radiation and compression is therefore 
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While the luminosity due to cooling by radiative losses 
diverges as r — > during inflow in the thermal peak, the 
radiative term in equation (3) nevertheless becomes pro- 
gressively smaller relative to the term for compressional 
heating which diverges even faster. As £ — > the gas 
temperature is altered only by compression, explaining 
the adiabatic nature of the flow inside the central ther- 
mal peak. 

3.1. Cooling Site Flows 

Cooling site solutions in which M — > are a well 
defined limit of radiating Bondi flows. In this limit 
a — > — 00 while M and A remain finite. Flows for large \a\ 
in Figure 1 remain essentially unchanged for \a\ > 100, 
so the \a\ = 100 profiles can be regarded as essentially 
equivalent to the (non-gravitating) cooling site flow. In 
cooling site flow the gas cools completely at some finite 
radius £ coo ; < 1 where all its entropy a = t/S 2 ^ 3 is ra- 
diated away. In these flows the total radiated emission 
within the sonic radius 



L c log(l/£ cooi ) 



is finite. 



4. CONFRONTING OBSERVATIONS 
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Of interest are potential sources of error when accre- 
tion rates are estimated by calibration with radiating and 
non-radiating Bondi flows. For this purpose it is neces- 
sary to convert from dimensionless to physical variables. 

4.1. Classical Bondi Flows 

In traditional Bondi accretion with 7 = 5/3 gas flows 
with constant entropy toward a mass point at the ori- 
gin. Very far from the mass point the gas is at rest with 
uniform density poc and temperature or sound speed c^, 
that define the adiabat for the entire flow. However, in 
Bondi solutions for a typical monatomic gas with 7 = 5/3 
the only flow from large radii to the origin must be fully 
subsonic, reaching the sound speed just at the origin, 
r s = 0. When 7 = 5/3 the energy integral for the steady 
flow is 



mal peak, having passed through a temperature mini- 
mum and a sonic point. The sonic radius, 
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If the gas velocity Uoo vanishes as r — > oo, the Bondi 
accretion rate is 



M = n(GM) 2 ^ 



(24) 



where the coefficient Poo/c^ is related to the inverse of 
the entropy factor, c^/poo oc (T/p 2 / 3 ) 3 / 2 . For a fixed 
central mass the accretion rate decreases with increasing 
gas entropy. When estimating the classical quadratic 
Bondi dependence, M oc M 2 , the entropy factor c^/poo 
in equation (24) is assumed to be independent of M. But 
it is easy to imagine astronomical correlations between 
entropy and M that are unrelated to the hydrodynamics: 
larger black hole masses are often found at the centers 
of more massive dark halos containing gas with higher 
entropy. 

The Bondi accretion rate M onto galaxy or cluster- 
centered black holes can be estimated from X-ray ob- 
servations of gas temperature and density profiles. Gas 
observed at some radius just larger than ~ GM/c s 2 , 
where the black hole potential begins to dominate, is 
assumed to be essentially at rest For example, Allen et 
al. (2006) estimated pa and Ta from (inward extrap- 
olated) observations at radius ta = GM/2ca 2 where 
ca = ("fkTA/ pnip) 1 / 2 is the local sound speed at ta 
and M is known from observations of the central galaxy. 
Since the entropy factor c s 3 / p is constant in Bondi flow, 
it can be observed at any radius. However, the Bondi 
flow morphology and equation (24) require Moo = 0. 

4.2. Radiating Bondi Flows 

In radiating Bondi flows radiation losses occur at all 
radii so the flow at large distances from the gravitat- 
ing mass point necessarily has gradients unlike the tra- 
ditional Bondi flow. For £ > 100, d log rj/d log £ = -2 
and d log r/d log £ = 1 are excellent approximations to 
radiating solutions in Figure 1 when £ > £ m j n . Conse- 
quently, the sound speed and gas density at infinity are 
not defined nor is this limit relevant in comparing with 
observation. 

Radiating accretion flows are inherently non-adiabatic. 
However, for \a\ 5s 1.5 radiating supersonic Bondi flows 
ultimately become adiabatic as £ — > in the central ther- 
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is likely to be small, not directly observable. Remarkably, 
the morphology of all radiating Bondi flows is determined 
by a single dimensionless parameter a. The accretion rate 
is 

M = a ■ 47rm 2 * , 26) 
p 5 A v ' 

where c s = u s is the sound speed at the sonic radius. 

Alternatively, by using equation (19) the radiating ac- 
cretion rate can be written in terms of the sonic radius 
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While the quadratic mass dependence in this less useful 
representation M oc M 2 resembles the classical Bondi 
relation, its coefficient a(15 — 6a)/r s cannot be directly 
observed and has a complex dependence on the flow vari- 
ables. 

Flows having the same M/M and T s (and therefore a) 
have identical temperature and density profiles T(r) and 
p(r) and for these flows the linear relation M oc M holds 
exactly. However, the relation M oc aMc 2 from equa- 
tion (26) does not in general imply that M oc M p with 
p = 1, nor does it imply that M and M are necessarily 
related with a simple power law. While accretion flows 
expressed in dimensionless variables, as in Figure 1, de- 
pend on the single parameter a, the same flows translated 
into dimensional variables u(r), T(r) and p(r) depend on 
any two of the parameters a, M and T s , assuming that M 
is independently known from observations of the central 
galaxy. 

A more general relationship between M and M can be 
established by considering intersecting flows. In particu- 
lar, by varying M and T s it is possible to construct radi- 
ating flows having different central masses M that share 
the same temperature T and density n e = p/(1.17m p ) 
at the same radius r. Such intersecting flows have differ- 
ent a and different morphologies. For example, Figure 
2 shows three radiating Bondi flows that intersect with 
identical gas properties. Only the subsonic regions of 
these flows are plotted, extending from the sonic radius 
r s at the left to 1 kpc at the right. The solid line profile 
is a reference solution with parameters M = 3 x 10 9 Mq, 
M = 0.3 M Q yr" 1 , and T s = 3 x 10 6 . K. The dashed 
lines are two radiating Bondi flows with lower central 
mass M/3 that intersect the reference flow at points 1 
and 2 where they have the same T and n e at r\ and r 2 . 
The upper dashed flow profiles (M/3)i in both panels 
of Figure 2 intersect at r\ — 152.2 pc (point 1) where 
T = 2.21 x 10 6 K and n e = 0.182 cm" 3 . The parameters 
of this flow are M = 10 9 M Q , M = 0.07917 M yr" 1 , 
and T s = 1.76 x 10 7 K. The lower dashed flow profiles 
(M/3)2 in Figure 2, based on parameters M = 10 9 M Q , 
M = 0.17929 Mq yr-\ and T s = 9.00 x 10 5 K, intersect 
at r = 624.3 pc (point 2) where T = 1.432 x 10 6 K and 
n e = 0.03568 cuT 3 . 
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Of interest are the differing values of M and M for the 
two flows involved at each of these intersections and in 
particular the mean power law variation M oc M p that 
approximately describes how the accretion rate varies 
with central mass for flows with central masses in the 
range M/3 — M, all of which intersect with the reference 
solution at the same point. From the two very specific 
intersections shown in Figure 2 at points 1 and 2 we find 
that M oc M p with p = 1.21 and 0.47 respectively. The 
exponent p differs from unity because the three intersect- 
ing flows in Figure 2 have different overall morphologies 
and values of a. The three flows in the upper panel of 
Figure 2 can be ordered - (M/3)i, M, (M/3) 2 - with 
decreasing r min /r s where r min is the radius of minimum 
gas temperature and the sonic radius r s is defined by 
the leftmost extent of each flow. A qualitatively similar 
morphological sequence for decreasing r min /r s = £ min is 
apparent in the central panel of Figure 1 for flows with 
\a\ = 0.01, 0.1 and 1. This is consistent with similarly in- 
creasing values of \a\ = 0.0137, 0.1008 and 0.6027 respec- 
tively that characterize the intersecting flows (M/3)i, M, 
and (M/3) 2 in Figure 2. 

The dotted lines (2M/3)i in Figure 2 show a third 
root through point 1 with parameters: 2M/3 = 2 x 10 9 
M Q , M = 0.1725 Mq yr" 1 , T s = 9.3 x 10 6 K and \a\ = 
0.02806. The variation of p between the three roots is not 
constant: p = 1.12 for (M/3)i ->■ (2M/3)i and p = 1.36 
for (2M/3)i — 5> M. Evidently, there is no general power 
law relation M oc M p that applies to all intersecting 
flows of this type, although local values of p are not far 
from unity. 

Also clear from Figure 2 is that the mass accretion rate 
M cannot be determined by observations of T and n e at a 
single radius. The parameters of radiating flows, M and 
T s , must be found with fits to a spatially extended set of 
observations that extend close to or within the radius of 
minimum temperature. The best-fitting flow parameters 
M and T s can be found with iterative integrations out 
from the sonic radius as we have done here. 

Three additional points: (1) Although we consider in- 
tersecting accretion flows having different M to explore 
the M oc M p relation, it must be stressed that not all 
flow solutions intersect nor do they intersect as shown 
in Figure 2. For example, subsonic entropy factor pro- 
files S(r) = T(r)/n e 2 / 3 (r) do not intersect when only the 
sonic temperature T s is varied. Subsonic entropy profiles 
in which only M varies typically intersect but, unlike 
Figure 2, the values of T and n e are not identical at 
intersections having the same entropy T/n e 2 / 3 . The un- 
usual flow intersections in Figure 2 with identical T and 
n e but varying mass M , T s and M reveal a relationship 
between a — a(T S} M /M) and the exponent p. Unlike 
traditional Bondi flows, the M(M) relation for radiating 
Bondi flows is not defined by the entropy T/n e 2 / 3 at a 
single radius of observation. (2) When M is known in ad- 
vance, the morphology parameter depends on the sonic 
temperature and accretion rate, a = a{T Sl M). As T s 
and M are varied to match observations of a single flow 
or when a group of similar flows are compared, the range 
of T s and M for successful fits is limited and correspond 
to a relatively small variation in a. When M and M are 



compared among these flows, we expect (3) We 

draw attention to the relative insensitivity of the index p 
to the morphology parameter a in the intersecting exam- 
ples illustrated in Figure 2. As \a\ increases from 0.0137 
to 0.6027, a factor of 44, p changes only from 1.21 to 0.47. 
In addition, for massive galaxy-centered black holes a is 
constrained by observations: \a\ ^ 3 is required to avoid 
total off-center cooling (eqn 28) which is not generally ob- 
served and \a\ cannot be very much less than 0.01 since 
temperature minima are rarely observed at large radii. 
This also suggests that p lies in a restricted region near 
unity. 

In Figure 3 we compare the accretion rate for the ref- 
erence solution M = 0.3 M© yr -1 with estimates of M 
using the classical Bondi procedure (eqn 24) . The Bondi 
mass flow rates M are determined from flow local gas 
temperature and density of the reference solution. It is 
seen that the classical Bondi solution considerably un- 
derestimates M in the region plotted. The decreasing 
classical Bondi M with radius mirrors the monotonically 
increasing entropy in the radiating solution. 

Equation (24) for Bondi flow is used extensively in 
computations of cosmological black hole accretion (e.g. 
Springcl et al. 2005; Johansson et al. 2009; Dubois et al. 
2010; Kim et al. 2011). But successful black hole growth 
is only possible if the Bondi accretion rate is enhanced 
by a large ad hoc dimensionless coefficient, while such a 
factor a appears naturally in radiating Bondi flows (eqn 
26). Bondi accretion is also invoked to estimate accretion 

rates M from X-ray observations of hot gas surround- 
ing galaxy-centered and cluster-centered black holes (e.g. 
Allen et al. 2006; Balmaverde et al. 2008; Vattakunnel et 
al. 2010). However, as with the cosmological computa- 
tions, X-ray observations of sources having greater feed- 
back power require mass accretion rates considerably in 
excess of the Bondi rate (Raffcrty et al. 2006; Hardcastle 
et al. 2007; McNamara et al. 2011). In addition to the 
large mass accretion rates implied by radiative emission 
from luminous black holes in AGNs and quasars, a com- 
parable accretion rate is required to provide the huge me- 
chanical power expended by cluster-centered black holes 
in expanding the cluster gas to match observed gas mass 
fraction profiles (Mathews & Guo 2011). 

In radiating Bondi flows the temperature dependence 
in this outer region £ > £ m j„, r oc £, is a reasonable 
match to gas temperature profiles in the central region 
of most galaxy groups and clusters. Moreover, the den- 
sity in £ > £, m i n varies asymptotically as 5 oc £~ 1//2 so 
the entropy is not constant but varies rather significantly 
with radius, a = t/S 2 ^ 3 oc £ 4 ^ 3 , at large £ in radiating 
Bondi flows. However, hot gas observed in more isolated 
elliptical galaxies can have dT/dr < 0. This can be un- 
derstood if the gas is confined by a halo of lower mass 
(Humphrey et al. 2006) or, at smaller radii, by compres- 
sional heating near the central black hole (Humphrey et 
al. 2008). 

Positive temperature gradients typically observed on 
kpc scales in galaxy groups and clusters are evidently 
imposed by entropy-losing cooling flows in the gravita- 
tional potential of the stars and the group/cluster dark 
matter on larger multi-kpc scales. These flows on larger 
scales have local enclosed masses greatly exceeding that 
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of the central black hole and contain gas with a larger 
virial temperature. Consequently, the gas flow rate M 
toward cluster-centered black holes depends not just on 
the black hole mass but on the huge mass of baryonic gas 
bound to the dark halo of the surrounding galaxy group 
or cluster. These large scale temperature and density 
gradients further confound estimates of the black hole 
accretion rate using the traditional Bondi recipe (eqn 
24). Once the innermost resolvable (or extrapolated) 
gas density and temperature are interpreted as and 
Too, the Bondi formula provides a unique local value of 
M = 7r(GM) 2 p 00 /c 02 for each M. By contrast, the corre- 
sponding relation for radiating Bondi flows, M cx aMT s , 
contains the additional dimensionless parameter a that 
can be adjusted to fit much larger mass flows M that ap- 
proach the central black hole from radiating gas bound 
to the surrounding group or cluster halos. 

Large values of the parameter \a\ can be regarded as 
a measure of large accretion rates that may cool signifi- 
cantly before reaching the central black hole. According 
to Figure 1, off-center cooling requires \a\ > 3 which, 
with equation (26), corresponds to larger accretion rates 

.Ml > 30 ( * ^ 1 J^^L_ M° 

1 '~ V!0~ 23 erg cm 3 s-V 10 7 K 1O 9 M yr ' 

(28) 

This is roughly consistent with observations of significant 
central star-formation in many X-ray luminous cool-core 
clusters (e.g. Donahue et al. 2010), although this oc- 
curs on scales larger than the sonic radius for the black 
hole. This inequality also expresses a biomdality between 
hot (|o| <> 3; \M\ <> 30 M Q yr^ 1 ) and cold (\a\ > 3; 
\M\ > 30) modes of accretion by central black holes. 
Radiating accretion flows that cool before reaching the 
central mass are those with small central mass M, lower 
temperature, or, as in the example above, large accretion 
rates M (Sarazin & White 1987; Quataert & Narayan 
2000). 

5. CONCLUSIONS 

Gas flowing near accreting black holes can only be 
observed with radiation that necessarily carries entropy 



away from the flow, in violation of the adiabatic assump- 
tion inherent in traditional Bondi flows. However, when 
radiative cooling is included, the accretion flow solu- 
tions are qualitatively different, not merely perturbations 
or limiting cases of the original Bondi flow. Neverthe- 
less, radiating solutions have widely varying morpholo- 
gies that depend on a single dimensionless parameter and 
they all must pass through a sonic point. Although the 
gas density (entropy) decreases (increases) monotonically 
with radius, many solutions exhibit a broad temperature 
minimum. Far from the mass point the temperature and 
entropy increase, allowing fits to be made to gas flow- 
ing subsonically in from a deeper potential well having 
a larger virial temperature. Some radiating accretion 
flows do not pass through a temperature minimum, but 
instead cool completely before reaching the central mass. 
In "cooling site" flows, in which there is no central mass, 
total cooling occurs at a non-zero radius within the sonic 
point. 

Of particular interest is the failure of traditional Bondi 
flows to accurately relate the mass accretion rate to the 
central mass. The traditional Bondi result M oc M p 
with p — 2 is certainly incorrect for radiating flows. The 
examples we consider indicate much smaller values p ~ 1 
in which the precise value of p is somewhat sensitive to 
the flow morphology. 

The idealized radiating Bondi flow we describe here, 
while certainly an improvement over the classic adiabatic 
Bondi flow, cannot be valid within some radius. Any 
small rotation in the initial gas must ultimately flow into 
a disk. Nevertheless, these idealized spherical radiating 
solutions may provide a useful framework for understand- 
ing more realistic flows that include rotation, magnetic 
fields, viscosity, feedback, etc. They may also be useful 
in evaluating the performance of central or non-central 
cooling in numerical computations. 

Our colleague Fabrizio Brighenti provided helpful ad- 
vice. We also acknowledge with thanks insights provided 
by the referee who challenged us to provide explicit ex- 
amples of non-quadratic relations between M and M. 
Studies of feedback gasdynamics in hot intracluster gas 
at UC Santa Cruz are supported by NSF and NASA 
grants for which we are very grateful. 
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Fig. 1. — Plot of radiating Bondi flows in dimcnsionlcss variables. Top: Dimcnsionlcss (square of the) velocity r;(£). Center: Dimcnsionlcss 
temperature t(£). Bottom: Dimcnsionlcss density <5(§). Where possible, each curve is labeled with its dimcnsionlcss parameter \a\. Profiles 
that are too crowded to label individually arc identified with the range in \a\ in order of decreasing variable n, r or <5. 
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Fig. 2. — Upper and Lower panels show respectively gas temperature and density profiles of the subsonic regions of three radiating 
accretion flows, each extending from 1 kpc at the right to the sonic radius r s at the left. For better visibility the gas density is multiplied 
by the radius in the lower panel. Solid lines show a reference solution for which M = 3 X 10 9 M Q , M = 0.3 M yr" 1 , and T s = 3 X 10 6 . 
K. The dashed lines show two additional flows (M/3)i and (M/3)2 for which the black hole mass M = 10 9 Mq is three times smaller. 
Parameters M and T B for the two dashed profiles have been chosen so that the intersections with the reference flow at points 1 and 2 at 
radius r have the same temperature and density as the reference solution. The dotted lines (2M/3)i show a third root through point 1 for 
an intermediate mass. The objective is to explore how M varies with M for flows observed having identical T and n e at some radius. 
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Fig. 3. — The solid line shows the constant mass accretion profile M = 0.3 for the reference radiating Bondi solution for M = 3 X 10 9 
Mq, M = 0.3 Mq, and T s = 3 X 10 6 . K. The dashed line shows the mass accretion rate for classic Bondi flows estimated with 
M = 7r(GM) 2 p 00 /c^ with the usual assumption poo = p(r) and Coo = c s (r). 



